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EFFECTS OF FLEXIBILITY ON LIFT AND 
PITCHING- MOMENT CHARACTERISTICS OF A 
SERIES OF LOW- ASPECT-RATIO WING-BODY 
COMBINATIONS AT TRANSONIC SPEEDS* 

By Robert V. Doggett, Jr., and A. Gerald Rainey 
Langley Research Center 

SUMMARY 

Measurements of lift and pitching moment were made on a series of low-aspect- 
ratio wing-body combinations of different stiffnesses in the Mach number range from 0.70 
to 1.10 and at angles of attack up to 30°. The configurations consisted of three different 
delta wings in combination with a single conical-cylindrical body. The three delta wings 
differed in leading- edge- sweep angle, aspect ratio, and area. Both rigid and flexible 
models were tested at Reynolds numbers per meter of 10.5 x 10® and 5.9 x 10®. The 
experimental results indicated that there is no appreciable effect of vehicle flexibility on 
the lift characteristics for the models tested. The pitching- moment characteristics were 
affected by flexibility; in general, increasing flexibility produced an appreciable deteriora- 
tion of the static longitudinal stability characteristics. Some of the experimental results 
were compared with some of the linear and nonlinear theoretical methods available. It 
is to be noted that in the theoretical calculations, the models were assumed to consist only 
of a triangular wing, effects of the fuselage being ignored. The linear theory showed good 
agreement with the experimental results in values of the lift- curve slope at the zero-lift 
condition. One of the nonlinear theories showed fair agreement with the experimentally 
determined lift characteristics for all three configurations. None of the analytical 
methods used were completely adequate for the prediction of the lift and pitching- moment 
characteristics throughout the entire angle-of-attack range of the investigation. 

INTRODUCTION 

The proposed use of low- aspect- ratio lifting vehicles for long-range hypersonic 
flight and for reentry into the earth's atmosphere has created a need for aerodynamic 

^Supersedes declassified NASA TM X-343 entitled "Investigation of Some Effects of 
Flexibility on the Lift and Pitching- Moment Characteristics of a Series of Low- Aspect- 
Ratio Wing-Body Combinations at Transonic Speeds." 



information to describe fully the behavior of such vehicles. One area where these vehi- 
cles have not received adequate attention consists of the aeroelastic problems of flutter, 
divergence, and deterioration of stability and control characteristics associated with 
elastic deformations. A preliminary study of these aeroelastic characteristics has indi- 
cated that this deterioration of stability and control characteristics might be very signif- 
icant. Most proposed ascent trajectories for these vehicles produce relatively high 
dynamic pressures in the transonic speed range. Since aeroelastic effects tend to become 
more pronounced at conditions where the product of dynamic pressure and lift-curve slope 
is a maximum, it appeared desirable to study the lift and pitching- moment characteristics 
of a series of low-aspect-ratio wing-body configurations in this often troublesome speed 
range. 

Accordingly, a program has been completed in the Langley 2 -foot transonic aero- 
elasticity tunnel where the lift and pitching moment have been measured on a series of 
full- span triangular-planform wing-body configurations varying in aspect ratio and stiff- 
ness. The measurements covered the Mach number range from 0.70 to 1.10 at angles of 
attack as high as 30°. The models were restrained at the trailing edge of the body by a 
balance which was attached to a conventional sting support system. The balance was 
shielded from the airstream by a fairing with a conical nose and cylindrical afterbody 
simulating a booster configuration. Thus, the lift and pitching- moment coefficients mea- 
sured may be interpreted directly as representing the aerodynamic loads imposed on a 
booster by such a forward- mounted vehicle. Also of interest, however, are the implica- 
tions contained in the data regarding the aeroelastic effects on the stability and control 
characteristics of free-flying vehicles of this type. The measured characteristics have 
been compared with those calculated by means of a relatively simple aeroelastic analysis. 

SYMBOLS 


A wing aspect ratio 

A ij deflection influence coefficient (deflection at ith point due to unit load at 

jth point) 

c wing mean aerodynamic chord 

c r wing root chord 

Ct lift coefficient, 

L ’ qS 

Cl lift- curve slope 
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C m pitching- moment coefficient, Pitching mome nt 

qSc 

C constant (see appendix) 

[d] differentiating matrix 

F aerodynamic force 

Fj aerodynamic force at jth point 

i,j integers 

M Mach number 

q dynamic pressure 

s = x tan e 

s 0 wing semispan 

S wing planform area 

x,y,z distances along X-, Y-, and Z-axes 

Ax element of width having its center at x = xj 

x_ n distance from leading edge of mean aerodynamic chord to center of pressure 

measured in fraction of mean aerodynamic chord, positive rearward 

Xj distance from origin to ith point 

Xj distance from origin to jth point 

z i deflection at ith point due to flexibility 

Zj deflection at jth point due to flexibility 

z'j deflection at jth point due to angle of attack 
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I 


a 


a. 


a e 


a i 


"ij 


angle of attack (measured at trailing edge of model) 


calculated angle of attack at station 1 


measured angle of attack at station 1 


rigid-body angle of attack at jth point 


calculated deflection at station 1 


measured deflection at station 1 


wing semiapex angle 


lambda function, Ay = 0 when i * j and Ay = 1 when i = j 


i] 


Matrix notations: 
square 


[] 

{ ) 


column 


APPARATUS AND TESTS 
Wind Tunnel 

The Langley 2-foot transonic aeroelasticity tunnel was used in this investigation. 
This tunnel is a slotted-throat single- return v/ind tunnel equipped to use either air or 
freon as the test medium at pressures from 1 atmosphere down to about 1/25 atmosphere. 
The tunnel is of the continuous-operation type and is powered by a motor-driven fan. 

Both test-section Mach number and density are continuously controllable. The present 
tests were made with freon as the test medium. Some of the characteristics of freon as 
a wind-tunnel test medium are discussed in reference 1. 


Models 

Model description .- Three series of models were tested. Details of the geometry 
of the models are shown in figure 1. The configurations consisted of three different low- 
aspect-ratio full-span delta wings in combination with a single conical- cylindrical body. 
Differences in the three series of wings were in leading- edge- sweep angle, aspect ratio, 
and area. The first series (hereinafter referred to as series A) had a sweep angle of 
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78.03° and an aspect ratio of 0.848. The second series (hereinafter referred to as 
series B) had a sweep angle of 83.95° and an aspect ratio of 0.424. The third series 
(hereinafter referred to as series C) had a sweep angle of 78.03° and an aspect ratio of 
0.848; however, the wing area was one-half that of the models of series A. The differ- 
ences between the models in any one series were in longitudinal bending stiffness. The 
individual models are hereinafter designated by two letters, and, in some cases, a number 
will be added. The first letter refers to the wing series; the second letter differentiates 
between the rigid models (the letter "R" is used) and the flexible models (the letter "F" 
is used). There were two flexible models of the A series having different levels of stiff- 
ness which are distinguished by the designations AF-1 and AF-2. Model AF-2 is the 
more flexible model. 

Model construction .- The models were constructed of aluminum, plastic, and a 
flexible plastic foam. The wings were molded from plastic. For the rigid models, the 
wings were impregnated with fiberglass to increase the wing stiffness. The fuselage was 
molded from a lightweight flexible plastic foam. The fuselages for the flexible models 
were cut transversely at approximately 2. 54 -centimeter intervals to reduce the contribu- 
tion of the fuselage stiffness to the total model stiffness. The cuts were covered with thin 
rubber sheet to preserve the aerodynamic contour and to prevent leakage. The primary 
contribution to the longitudinal bending stiffness of the models was provided by an 
aluminum-alloy spar which was bonded to the root chord of the wing. Variations of the 
stiffness were obtained by variations of the dimensions of the spar. Two transition strips 
of No. 60 carborundum grains were applied near the nose of the rigid models to simulate 
the roughness of the first two rubber- covered cuts on the flexible models. A photograph 
of a typical flexible model of series A with part of the fuselage removed is shown in 
figure 2. 

Presented in table I is the measured deflection influence coefficient matrices for 
the flexible models. These coefficients are given as deflection in millimeters per newton 
for each of the 10 stations on the model shown in figure 3. Since the model was mounted 
on a strain-gage balance during the wind-tunnel tests, the influence coefficients were mea- 
sured with the models mounted on the balance. It should be noted that the balance was 
relatively stiff. The flexibility of the balance may be conveniently expressed as a rota- 
tional spring constant. This constant was determined to be about 3616 meter- 
newtons/radian with the effective axis of rotation located 5.24 centimeters rearward of 
the model trailing edge. As an example of the effect of the balance on the overall flexibil- 
ity of the model-balance system, the deflection at station 3 of model AF-1 due to a load 
at that point was increased by about 10 percent because of the balance flexibility. It 
should be noted that the rigid models were not infinitely stiff and perhaps would be more 
properly described as being very stiff in comparison with the flexible models. 
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Model support s ystem.- The models were cantilever-mounted on a three- component 
strain-gage balance which was in turn attached to a support sting. A line drawing showing 
the planform of the model support system is shown in figure 4. The balance was shielded 
from the airstream by a conical- cylindrical fairing, the geometry of which is representa- 
tive of the forward portion of a typical rocket booster system which would normally be 
attached to a full-scale vehicle. A line drawing of the fairing is shown in figure 5. A gap 
of approximately 0.159 centimeter was left between the trailing edge of the model and the 
balance shield. The sting support was attached to a circular turntable which formed part 
of the tunnel wall. By rotating the turntable, the model angle of attack was varied. A 
photograph of a typical model mounted in the test section is shown in figure 6. 

INSTRUMENTATION AND DATA REDUCTION 

The forces and moments acting on the models were measured by means of a three- 
component strain-gage balance. The normal-force, chord-force, and pitching- moment 
components of the balance were designed for maximum measurable loads of 534 newtons, 
356 newtons, and 27 meter- newtons, respectively, with an accuracy of ±1/2 percent of the 
maximum load. 

At relatively small angles of attack the measured chord forces were small and, 
within the precision of the balances used, could not be determined with sufficient accuracy 
to justify their presentation. The pitching- moment measurements appeared to be affected 
by changes in the temperature of the balance. These effects caused some drift in the zero 
reading of the balance; however, since the calibration of the balance was essentially unaf- 
fected by temperature, the slopes or shapes of the pitching- moment curves were correct 
because all runs were made at essentially constant temperature. 

All the force and moment data have been reduced to coefficient form. The ref- 
erence length and reference area used were the wing mean aerodynamic chord and the 
total wing area, respectively. The pitching- moment coefficient is referred to an axis on 
the surface of the wing, parallel to the wing trailing edge, and located at 42 percent of the 
mean aerodynamic chord. (See fig. 1.) The angle of attack a is defined as the angle at 
the trailing edge. 

TEST CONDITIONS 

The aerodynamic lift forces, drag forces, and pitching moments were determined 
at Mach numbers of 0.70, 0.80, 0.85, 0.90, 0.95, 1.00, 1.05, and 1.10 and at angles of 
attack from -4° to 30°. All the models were tested at a Reynolds number per meter of 
approximately 10.5 x 10®. Tests were also made on the models of series A at a Reynolds 
number per meter of approximately 5.9 x 10®. Shown in figure 7 is the variation of the 
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test-section dynamic pressure with Mach number for the two test Reynolds numbers. 
Since there were slight variations in the tunnel stagnation pressure between tests, the 
dynamic pressure at a given Mach number varied from run to run. The data in the fig- 
ure give the maximum range covered for all runs. This variation in the tunnel stagna- 
tion pressure caused a maximum deviation in the Reynolds numbers of approximately 
2 percent. 


EXPERIMENTAL RESULTS AND DISCUSSION 
Lift Data 

Rigid models .- Presented in figure 8 is the variation of the lift coefficient with 
angle of attack for all the rigid models tested. In order to facilitate presentation of the 
data, staggered scales have been used in many of the figures and care should be taken in 
identifying the zero axis for each curve. In figure 8(a), data from tests on model AR at 
Reynolds numbers per meter of approximately 10.5 X 10® and 5.9 x 10® are presented. 
There appears to be no appreciable effect of Reynolds number on the variation of the lift 
coefficient with angle of attack for the rigid model of series A. Slight differences are 
noted in the data taken at the two Reynolds numbers, but these differences are within the 
experimental error. Consequently, in subsequent comparisons of lift data for the flexible 
models, the effects of Reynolds number are assumed to be negligible. Presented in fig- 
ure 9 is a comparison of the lift- coefficient data for the rigid models of the three dif- 
ferent configurations at several Mach numbers (M = 0.70, 0.90, 1.00, and 1.10). A gen- 
eral comparison of the data in figure 9 indicates that all the models show essentially the 
same type of variation of the lift coefficient with angle of attack. 

Flexible models .- Presented in figure 10 is the variation of the lift coefficient with 
angle of attack for the flexible models. Also included in the figure are the faired curves 
from figure 8 for the rigid models. On comparing the data for the flexible models with 
that for the rigid models, it is seen that the lift data are, for all practical purposes, 
unaffected by variations in model flexibility. It should be noted that the angle of attack 
is measured at the trailing edge. 


Pitching- Moment Data 

Rigid models .- Presented in figure 11 is the variation of the pitching- moment coef- 
ficient with lift coefficient for the rigid models. As is seen from figure 11(a), there are 
some differences in the data taken at the two Reynolds numbers for model AR. These 
differences may be associated with viscous effects; however, some of the difference may 
be attributed to inaccuracies in the experimental method or to the fact that the rigid 
model was not infinitely stiff. Since the tests at two Reynolds numbers involve two sets 
of values of dynamic pressure (shown in fig. 7), effects of deformation on the data may 
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have been present. The general trend shown by the pitching- moment data for positive 
values of the lift coefficient is an almost linear decrease until some value of the lift coef- 
ficient, depending on the Mach number, is reached and then the pitching- moment data show 
an increase in value. The slope of the linear portion of this variation tends to become 
more negative with increasing Mach number. 

Flexible mode ls.- Presented in figure 12 is the variation of the pitching- moment 
coefficient with lift coefficient for the flexible models. An examination of these data 
shows that the trends are similar to those described for the rigid model. 

A more direct indication of the effects of flexibility on the pitching-moment charac- 
teristics is given in figure 13 where the variation of pitching- moment coefficient with lift 
coefficient is shown for the models of the A series having three different stiffnesses. 

When the data are examined from the standpoint of the stability and control characteris- 
tics of a free-flying vehicle, several interesting aspects are noted. The rigid model would 
be stable about the chosen center-of-gravity axis at 42 percent of the mean aerodynamic 
chord for all Mach numbers. The rigid model indicates a pitch-up instability at high lift 
coefficients for the lowest Mach number (M = 0.70). At positive lift coefficients and at 
the lower Mach numbers, the models of increasing flexibility show a systematic decrease 
in the lift coefficient at which pitch-up occurs. In addition, the most flexible model shows 
an unstable pitching- moment characteristic at negative lift coefficients for all Mach num- 
bers. These rather serious effects of flexibility on the static longitudinal stability char- 
acteristics indicate that the aeroelastic characteristics of highly flexible vehicles of this 
type should be given careful consideration. 

COMPARISONS WITH THEORY 
Lift Data 

Several theories (refs. 2 to 8) are available for the prediction of the lift for low- 
aspect-ratio triangular wings. The theories of references 2 to 5 are developed for a rigid 
wing in a steady incompressible flow. In reference 6 a linear aerodynamic theory has 
been developed for a wing which can deform elastically in the camber direction. The 
theoretical method presented in reference 7 is similar to that of reference 6 but is more 
general since both steady and unsteady aerodynamic forces are considered. The subsonic 
lifting- surface theory presented in reference 8 has application to both the steady and 
unsteady case and also permits effects of model deformations on the aerodynamic forces 
to be taken into account. 

The theory of reference 2 is based on the idealization of two-dimensional incom- 
pressible flow and gives for the lift coefficient at small angles of attack 
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( 1 ) 


C L = 


irAa 

2 


References 3 to 5 are extensions of the work of reference 2 to include nonlinear effects 
of viscosity. In these theories the expressions for the lift-force coefficient consist of a 
linear term (the result given by ref. 2) plus a nonlinear viscous term. However, by 
making slightly different assumptions as to the nature of the flow field in the vicinity of 
the wing, different expressions for the viscous term were determined. The expressions 
for the lift coefficients given by these references are as follows: 


From reference 3, 


from reference 4, 


and from reference 5, 


irAa. A(nat)3/2 
2 4 


C L 


^M+ttA VV / 3 


C L 


nA a 

2 


+ ^/V / 3 



( 2 ) 

(3) 

(4) 


It is to be noted that all calculations were made by considering the models to con- 
sist only of a triangular wing, effects of the fuselage being ignored. A comparison of the 
experimentally determined lift coefficients for the rigid models with the theoretical 
values obtained from equations (1) to (4) is presented in figure 14. Only experimental 
data for M = 0.70 are included in the figure since these data are typical of all data 
obtained. Also included in figure 14(a) are the results obtained by using a subsonic 
lifting- surface treatment. This method is developed in reference 8 for the unsteady case 
by using a spanwise numerical integration of the kernel function of unsteady, three- 
dimensional compressible flow. By considering the case where the reduced frequency 
equals zero, results for steady flow are obtained. The curve presented has been calcu- 
lated for M = 0.70. As is seen from figure 14(a), none of the theories are satisfactory 
for predicting the variation of the lift coefficient with angle of attack throughout the entire 
range of angle of attack for the models of series A and C. At low angles of attack, the 
linear theory and the lifting- surface theory (refs. 2 and 8) show the best agreement with 
the experimental data. At higher angles of attack ( a > 10°), the theory of reference 3 
gives the best prediction for model AR and the lifting- surface theory (ref. 8) gives the 
best prediction for model CR. In figure 14(b) the theory of reference 3 satisfactorily 
predicts the variation of the lift coefficient with angle of attack for model BR throughout 
the entire angle-of-attack range. 
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Presented in figure 1 5 is the variation with Mach number of the experimental lift- 
curve slope at the zero- lift condition for the models of all three series. The experimen- 
tal values were determined by numerically differentiating the experimental lift data with 
a five-point differentiating scheme. Also included in the figure are the corresponding 
theoretical values obtained from the theory of reference 2. Included in figure 15(a) is 
the calculated lift- curve slope obtained by the method of reference 8. As is seen from 
the figure, the results obtained from the linear low -aspect- ratio aerodynamic theory are 
in good agreement with the experimental data. The three-dimensional theory of refer- 
ence 8 predicts a value of the lift- curve slope somewhat higher than the value found 
experimentally . 

It is observed that the usual relatively large deviations in lift- curve slope found for 
most configurations in the transonic speed range were not encountered in these tests. 

This result is believed to be due to the slenderness of the configurations tested. It is of 
interest that the theory of reference 8, which takes account of Mach number, also shows 
the same trend as the experiments. 

In reference 6 linear low-aspect-ratio aerodynamic theory has been applied to a 
low-aspect-ratio triangular wing which is allowed to deform elastically in the camber 
direction. The total lift force was found to be independent of the deformations and depen- 
dent only on the slope at the trailing edge. Thus, the lift on the flexible wing at a given 
angle of attack, measured at its trailing edge, would be the same as the lift on a rigid 
wing. Apparently, the increased lifts generated by the higher local angles near the nose 
are compensated for by the negative lift forces associated with the curvature on the wing. 
Although linear theory is inadequate for the prediction of the magnitudes of the lift coeffi- 
cient throughout the test angle- of- attack range, the experimental results are consistent 
with the results which would be expected from the method of reference 8 in that no appre- 
ciable effect of flexibility was found in the lift data. (See figs. 10 and 15.) 

Center-of-Pressure Data 

The location of the aerodynamic center of pressure as a fraction of mean aerody- 
namic chord for model AR is presented in figure 16. Also included is the center-of- 
pressure location as predicted by the theories of references 2 and 8. Both theories pre- 
dict a more rearward location than was found experimentally. Both the experiment and 
the lifting- surface theory (ref. 8) show a rearward movement of the center of pressure 
with increasing Mach number. 


Pitching- Moment Data 

Figure 17 presents a comparison of some experimental and calculated variations of 
the pitching- moment coefficient with lift coefficient for the three series of models. 
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Calculations were made for all the configurations by using the theory of reference 2. 
Calculations were also made for one of the flexible models of series A by using the linear 
low-aspect- ratio aerodynamic theory and allowing the model to deform elastically in the 
camber direction (ref. 7). The deformed shape used in this calculation is presented in 
figure 18. As is seen from figure 17(a), the calculations for the rigid model based on 
reference 2 do not satisfactorily predict the variation of the pitching- moment coefficient 
with lift coefficient for any of the three configurations. Of course, this linear theory 
could not be expected to predict the pitch-up tendency at high lift coefficients. In addi- 
tion, some of the discrepancy between theory and experiment may be due to the omission 
of the effects of the body in the calculation. In figure 17(b) it can be seen that the calcu- 
lation for the flexible model shows much better agreement with experiment than did the 
calculations for the rigid models; however, the calculation for the flexible model cannot 
be considered to yield an adequate estimate of the pitching- moment characteristics. 

Elastic Deformations 

A comparison of calculated and measured normalized shapes for model AF-1 is 
presented in figure 18. The type of deformed shape obtained on the flexible models was 
basically a bending in the camber direction. Some significant bending in the spanwise 
direction was also observed at high angles of attack. The experimental deflection shapes 
were determined by taking a double- exposure photograph of model AF-1 which had been 
painted black with very thin stripes of white running in the camber direction. One expo- 
sure was made with the wind off and the model at 0° angle of attack; whereas, the second 
exposure was taken at the desired test condition. By measuring the deflections of the 
white lines, the deformed shape of the model was determined. The calculated deflection 
shapes were determined by using a matrix iteration technique employing the linear low- 
aspect-ratio aerodynamic theory of reference 7 and measured deflection influence coef- 
ficients. This analysis is developed in detail in the appendix. The shapes presented in 
figure 18 have been normalized to the value of the deflection at station 1. (See fig. 3 for 
station 1 location.) The actual deflection at station 1 is indicated in figure 18. The 
values of both the calculated and measured angles of attack at station 1 are also tabulated 
in figure 18. It is seen in figure 18 that the general shapes of the experimental and theo- 
retical curves compare favorably. However, the theory predicts a higher total deflection 
than was found experimentally. 

As a matter of reference, the divergence dynamic pressures were calculated for all 
the flexible models tested by using a technique similar to that developed in reference 9. 
This technique is developed in the appendix. These values of the dynamic pressure were 
135 kN/n^ 71.6 kN/m2, and 221 kN/m2 for models AF-1, AF-2, and BF, respectively. 

At the maximum test conditions, models AF-1, AF-2, and BF were tested at approximately 
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9 percent, 17 percent, and 6 percent of their calculated divergence dynamic pressures, 
respectively. 


CONCLUDING REMARKS 

Measurements of lift and pitching moment were made on a series of low-aspect- 
ratio wing-body combinations of different stiffnesses in the Mach number range from 0.70 
to 1.10 and at angles of attack up to 30°. The configurations consisted of three different 
delta wings in combination with a single conical- cylindrical body. The three delta wings 
differed in leading- edge- sweep angle, aspect ratio, and area. Both rigid and flexible 
models were tested at Reynolds numbers per meter of 10.5 x 10® and 5.9 x 10®. 

The experimental results indicated that there is no appreciable effect of vehicle 
flexibility on the lift characteristics for the models tested. The pitching- moment char- 
acteristics were affected by flexibility; in general, increasing flexibility produced an 
appreciable deterioration of the static longitudinal stability characteristics. 

Some of the experimental results were compared with some of the linear and non- 
linear theoretical methods available. It is to be noted that in the theoretical calculations 
the models were assumed to consist only of a triangular wing, effects of the fuselage 
being ignored. The linear theory showed good agreement with the experimental results 
in values of the lift- curve slope at the zero-lift condition. One of the nonlinear theories 
showed fair agreement with the experimentally determined lift characteristics for all 
three configurations. None of the analytical methods used were completely adequate for 
the prediction of the lift and pitching- moment characteristics throughout the entire angle- 
of- attack range. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., July 12, 1960, 

124-06-05-04-23. 
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APPENDIX 


TECHNIQUE EMPLOYED IN CALCULATING DEFORMED SHAPES 

A method of analysis was developed for calculating the deformed shapes of the 
flexible models by using the aerodynamic forces obtained from the theory presented in 
reference 7 for low- aspect- ratio triangular wings and measured deflection influence coef- 
ficients. The model is represented structurally and geometrically as shown in the fol- 
lowing sketch: 


Z 



The wing in its neutral position is assumed to have its mean camber surface lying in the 
XY-plane of the X,Y,Z coordinate axes system with the wing apex at the origin of the 
coordinate system and the root chord coinciding with the X-axis. The Z-axis is taken as 
positive upward, and a wind of constant velocity and inclination a to the XY-plane ema- 
nates from the negative x-direction. The trailing edge of the wing is considered to be 
built in and the wing is allowed to only have deformations in the camber direction. 

For any type of loading, the elastic deflection at any point x = x^ on the wing is 
defined by 

z i = A il F l + A i2 F 2 + A i3 F 3 • • • + A ij F j ( A1 ) 

where A^j is an elastic deflection influence coefficient defined as the deflection at posi- 
tion x = Xj due to a unit load at position x = xj. For the present analysis the wing was 
divided into 10 chordwise segments. Generalizing equation (Al) to include all the control 
points on the wing leads to the following deflection influence coefficient equation: 

■ NO'S} < A2 > 
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For the case where a = 0°, the aerodynamic force per unit chord as given by ref- 
erence 7 is 

<1F _ f„ n 2,Y v 2 (a: 


By considering an element of width Ax having its center at x = xj, the total aerody- 
namic force on this element is given by 


Fj = 2„(A*)q tan 2 s Xj 2 (0} + 2x i(S)j 


Rewriting equation (A4) in matrix notation for the entire wing gives the following equation: 


{Fj} = 2ir(Ax)q + 2x j(f ) 


With the use of the expressions Ay = 1 when i = j and A^j = 0 when i * j, equa- 
tion (A5) can be rewritten in the form 

<?i> = 2,,<Ax)q tan2e (^‘J x i 2 ]((^j) + 2 M 


By determining a differentiating matrix [dJ such that 




equation (A6) may be written in terms of slope only. The resulting equation is 


0i} = 2 *< Ax > q + “Mil"®) 


Since all the premultipliers of the matrix i n equation (A8) with the exception 

of q are constants peculiar to the particular configuration, the aerodynamic force may 
be expressed as 


(Fj> - q[c] 


Again, if the differentiating matrix [ D ] is employed so that 


i i ■■■■■■ mmm i mb ■■■■■■ i laHlill nimiiiii II Ml 



APPENDIX 


(iS),}' BCD 

the aerodynamic force becomes 

{Fj}=q[c][D]{ Z J> 

On substituting equation (All) into the deflection influence coefficient equation (eq. (A2)), 
the final equation for the case a = 0° is obtained: 

(*i} = q [Alj][c][D]{ij} (A12) 


(A10) 

(All) 


By iterating equation (A12) for the dominant root, the divergence dynamic pressure is 
obtained. This result is essentially the same as that presented in reference 

For the case a * 0°, there is an aerodynamic force associated with the angle of 
attack in addition to the force due to the elastic deformations. The force for the wing 
acting as a rigid body at an angle of attack may be expressed in matrix notation by the 
following equation: 

{Fj} = 27r(Ax)q tan 2 e J 2 [>-ij x jjJ j} (A13) 

Since 

[ D ){“i} = 0 


equation (A13) may be rewritten in the form of equation (A8) 

(Fj) = 27r(Ax)q tan^jjxijX^jD] + 2(\ ij xj] {aj} 
With the use of the matrix [jcfj and the relationship 

' C D K z 'i) 


(A14) 


equation (A14) becomes 


{ F i> - q [ c ][ D ]{ z ’i} < A15 > 

By adding equations (All) and (A15) and substituting the result into the deflection influence 
coefficient equation (eq. (A2)), the final equation for the case a * 0 is obtained: 

{5 O’ = q [ A ii][ c ][ D ]{ z i + z 'i} ^ A16 ^ 

For the case a= 0° equation (A16) reduces to equation (A12). Equation (A16) is iter- 
ated for the deformed wing shape by using the desired initial angle of attack and dynamic 
pressure. 
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TABLE I.- DEFLECTION INFLUENCE COEFFICIENT MATRICES 


(a) Model AF-1 










00 


TABLE I.- DEFLECTION INFLUENCE COEFFICIENT MATRICES - Continued 


(b) Model AF-2 






Aij, 

mra/N, at j of 

- 




i 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

i 

1.1871 

0.8160 

0.5773 

0.4003 

0.2729 

0.1759 

0.0999 

0.0422 



2 

.8160 

.6446 

.4819 

.3403 

.2370 

.1530 

.0868 

.0377 

.0160 

.0086 

3 

i .5773 

.4819 

.3866 

.2861 

.2010 

.1302 

.0736 

.0331 

.0148 

.0080 

4 

.4003 

1 

.3432 

.2861 

.2290 

1 

.1650 

.1073 

.0605 

.0286 

.0137 

.0074 

5 

.2729 

.2370 

.2010 

! .1650 

.1290 

.0845 

.0474 

j .0240 

.0126 

.0068 

; 6 

.1759 

.1530 

.1302 

.1073 

.0845 

.0588 

.0343 

| .0194 

.0114 

.0063 

7 

1 .0999 

! .0868 

.0736 

.0605 

.0474 

i .0343 

| 

1 .0211 

i .0148 

.0103 

.0057 

8 

.0422 

.0377 

.0331 

.0286 

.0240 

j .0194 

.0148 

| .0103 

.0086 

.0051 

9 

.0171 

.0160 

.0148 

.0137 

.0126 

! .0114 

j .0103 

.0091 

.0080 

.0046 

i 10 

.0091 

1 

.0086 

, .0080 

.1 

.0074 

.0068 

.0063 

.0057 

; .oo5i 

.0046 

j .0040 
.1 











TABLE L- DEFLECTION INFLUENCE COEFFICIENT MATRICES - Concluded 


(c) Model BF 






A^j, mm/N, at j of 

- 




i 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

i 

0.7263 



0.2615 

0.1839 

0.1193 

0.0674 

0.0308 

0.0166 


2 

.5048 

.3963 


.2238 

.1582 

.1028 

.0594 

.0274 

.0154 

.0086 

3 

.3637 



.1861 

.1325 

.0862 

.0514 

.0240 



4 

.2615 

.2238 

.1861 

.1485 


.0697 

.0434 

.0206 


.0074 

5 

.1839 

.1582 

.1325 


.0811 

.0531 

.0354 



.0068 

6 

.1193 


.0862 

.0697 

.0531 

.0365 

.0274 


.0108 

.0063 

7 

.0674 


.0514 

.0434 

.0354 

.0274 

.0194 


mm 

.0057 

8 

.0308 


.0240 


.0171 

.0137 

.0103 



.0051 

9 

.0166 


.0143 

.0131 

1 





.0046 

10 

.0091 



.0074 


K3I 































•Tunnel wall 



Figure 4.- Planform of model support system. 
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(a) Model AR. 

Figure 8.- Variation of lift coefficient with angle of attack for rigid models. 
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(b) Model BR. Reynolds number per meter ~ 10.5 * IQ 6 . 


Figure 8.- Continued. 
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(a) Model AF-1. 

Figure 10.- Variation of lift coefficient with angle of attack for flexible models of series A and series B, and a comparison with 

rigid-model curves from figure 8. 
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(b) Model AF-2. 


12.- Continued. 
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(c) Model BF. Reynolds number per meter ~ 10.5 x 10 6 . 
Figure 12.- Concluded. 



Figure 13.- Comparison of variation of pitching-moment coefficient with lift coefficient for models of series A. 

Reynolds number per meter - 10.5 x 10 6 . 
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(b) Series B. 
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(c) Series C. 

Figure 15.- Variation of lift-curve slope at zero- lift condition with Mach number. Reynolds number per meter ~ 10.5 x 10 6 . 












(a) Rigid models. 

Figure 17.- Comparison of experimental and calculated variations of pitching-moment coefficients with lift coefficient. 
M = 0.70; Reynolds number per meter ~ 10.5 x 10 6 . 
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Figure 18.- Comparison of calculated and measured normalized deflection shapes for model AF-1. Reynolds number per meter ~ 10.5 x hA 
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